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ABSTRACT 

We fit the colour-magnitude diagrams of stars between the zero-age main-sequence and 
terminal- age main sequence in young clusters and associations. The ages we derive are 
a factor 1.5 to 2 longer than the commonly used ages for these regions, which are 
derived from the positions of pre-main-sequence stars in colour-magnitude diagrams. 
From an examination of the uncertainties in the main-sequence and pre-main-sequence 
models, we conclude that the longer age scale is probably the correct one, which implies 
we must revise upwards the commonly used ages for young clusters and associations. 
Such a revision would explain the discrepancy between the observational lifetimes of 
proto-planetary discs and theoretical calculations of the time to form planets. It would 
also explain the absence of clusters with ages between 5 and 30Myr. 

We use the r 2 statistic to fit the main-sequence data, but find that we must make 
significant modifications if we are to fit sequences which have vertical segments in the 
colour-magnitude diagram. We present this modification along with improvements to 
methods of calculating the goodness-of-fit statistic and parameter uncertainties. 

Software implementing the methods described in this paper is available from 
|http://www. astro. ex. ac.uk / people / tim n / tau-squared/| 

Key words: stars: formation - stars: pre-main-sequence - methods: statistical - open 
clusters and associations: general - stars: fundamental parameters - stars: early-type. 



1 INTRODUCTION 

The ability to determine the ages of pre-main-sequence 
(PMS) stars is crucial for advancing our understanding of 
the early phases of stellar evolution. There are two key ap- 
plications. Firstly, and perhaps most obviously, we need stel- 
lar ages if one is to carry our experiments such as tracing 
the evolution of stellar angular momentum, or following the 
fraction of stars with proto-planetary discs as a function 
of time. Secondly, for PMS stars, the conversion from ob- 
servables such as temperature and luminosity into mass is 
highly age dependent, making accurate ages vital for deter- 
mining the mass function. The primary method of deter- 
mining the ages required for these studies is to compare the 
observed properties of PMS stars with models. The most 
easily accessible observables are a star's temperature and 
luminosity, since they can be measured from its colours and 
magnitudes. The problem is that for the same colours and 
magnitudes different models can predict ages which differ 
by a factor two, and even the same models will predict dif- 
ferent ages depending on which colours and magnitudes are 
used. This makes meaningful comparisons between the ages 
quoted in the literature for clusters or associations at best 
difficult, and often impossible. It was these problems which 
led us to devise a model-independent age ordering of young 
clusters and associations based on their colour-magnitude 



diagrams ( Mayne et al.|[2007 ). For PMS stars the primary 
age diagnostic is based on the fact that stars fade as they 
get older and contract towards the main sequence (MS). We 
used this movement of the sequence towards progressively 
fainter magnitudes to derive an age ordering, although to 
do so we also had to measure a consistent set of distances, 
which we derived from the more massive stars which have 
already reached the MS ( |Mayne fe Naylor||2008| ) . 

Whilst an age ordering such as ours is useful, for ex- 
ample it has showed unambiguously that different clusters 
take different times to reach the same disc fraction or an- 
gular momentum distribution, for quantitative work an ab- 
solute scale is required. For PMS clusters and associations 
there are several usable age indicators, each of which re- 
lies on comparing stellar properties with models. For this 
reason it is best to group them according to the underly- 
ing physics. First is the contraction of PMS stars as they 
approach the MS. As pointed out above, and discussed at 
length in Mayne et alT1(|2Q07|) these "contraction" or "PMS" 



ages are highly model dependent, and given the current dis- 
agreements between the models cannot yield an absolute age 
scale. Although most stars in a young cluster or association 
are in the PMS phase, the evolution of the most massive 
stars proceeds so fast that they may not only have reached 
the MS, but evolved beyond it. This gives us access to two 
more age measures. First, having reached the MS, stars move 
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redwards and to higher luminosities away from the zero-age 
main sequence (ZAMS), due to the increasing helium con- 
tent of their cores. This movement continues until the point 
of core hydrogen exhaustion, when the star has reached the 
reached the terminal age MS (TAMS or MS turn-off). Fi- 
nally, after the turn-off, the post-main-sequence evolution 
is driven by the burning of heavier elements which leads to 
much more rapid movement in the CMD. This relatively 
high velocity in the CMD means that post-main-sequence 
evolution has the potential to give precise ages. However, 
for young galactic clusters the paucity of stars in this re- 
gion of the CMD means such an age can depend on just 
one star, and such ages are rightly treated with some scep- 
ticism. Conversely, the main-sequence evolution (from the 
ZAMS to the turn-off) has a larger number of stars, but the 
movement is often subtle, and using the normal technique 
of simply plotting isochrones over the data leads to large 
uncertainties in age, and to questions over objectivity. How- 
ever, we have been developing a method of making objective 
fits to colour-magnitude data, which should allow us to un- 
lock the information in this stage of a star's evolution. The 
technique, called r 2 fitting, can be viewed as an extension of 
X 2 to data points with uncertainties in two or more observ- 
ables, and to models which are distributions (not just lines) 
in the data space. 

The aim of this paper is to apply the r 2 fitting tech- 
nique to the main-sequence evolution of young stars, and 
use the resulting ages to create a revised age scale for PMS 
stars. Surprisingly, this leads to a significantly older ages 
than the commonly used contraction ages, a result which 
we will discuss in Section [TT] To derive this result we first 
have to update our statistical techniques originally described 
Naylor & Jeffries (2006), since, as we discuss in Section [T] 



the technique will not work for the isochrones we wish to fit. 
We therefore lay out the changes which need to be made by 
following an example through fitting (Section^, testing the 
goodness of fit (Section [6J and determining the uncertain- 
ties in the derived parameters (Section [7]). Before doing so, 
however, we discuss the data and models we use (Sections 
[2] and Sections |3| . We deal with the effects of interstellar 
extinction in Section [8] and the details of each cluster in 
Section |9] We draw all the results together in our discussion 
in Section [TT] 



2 THE DATA 

To compare a set of ages derived from MS evolution with 
contraction ages we need a sample of clusters and associ- 
ations which have contraction ages, and for each of which 
data are available for MS fitting. Our sample is, therefore, 
based on the groups we placed in age order using the PMS in 
|Mayne &; Naylor| ( |2008| . Clearly, for each of these groups we 
require stars in the appropriate mass range to show signifi- 
cant MS evolution, but we also require extinctions and reli- 
able distance measurements. UBV photometry can provide 
all three of these. First the U — B / B — V diagram provides 
extinctions. Second, the upper part of V/B — V diagram is 
age sensitive, tracing the evolution of stars from the ZAMS 
to the turn-off. Finally, in the age range of interest the lower 
mass stars are still close to the ZAMS, and the sequence 
turns redwards, making it ideal as a distance measure. Fur- 



thermore, the UBV photo-electric system is very consistent 
and well characterised. However, to ensure we maintain the 
highest level of consistency we have restricted ourselves as 
far as possible to the data of Johnson and collaborators, pri- 
marily taken in the 1950s and 1960s. As we shall show later, 
the qualit y of these data whe n combined with the transfor- 
mations oflBessell et al. 1(1 19981 is impressive, giving r values 



which mean the model is a good fit to the data. Clearly we 
wish to avoid PMS stars contaminating our sample at faint 
magnitudes and red colours, and so for most objects we ap- 
ply a cut in observed B — V which roughly corresponds to 
(B-V)o <0.0. 

Most of the datasets we use have robust uncertainties 
derived from comparisons of many measurements of stars. 
This presents us with a problem, as the quoted uncertainties 
in colour are always smaller than those in magnitude. Con- 
ventional error analysis yields a correlation between, say, V 
and B — V, and in previous work we have always been care- 
ful to include that correlation when modeling the uncertain- 
ties. The starting point for such an analysis is that V and 
B are measured independently, and so the uncertainties in 
V and B — V are SV and ^/SV 2 + SB 2 respectively. Such an 
analysis also leads to the conclusion that the uncertainty in 
B — V must be larger than that in V, in direct contradiction 
to the quoted uncertainties for most of the data presented 
here. This is because it is not photon statistics which are 
the driver of the uncertainties, but changes in the trans- 
parency. In this work, we therefore model the uncertainties 
as uncorrelated. 



3 THE MODELS 

Although we will try other models later, we begin by using 
"Geneva-Bessell" isochrones. For the stellar interior we fol- 



low the suggestion of Lejeune & Schaerer (2001), and use 
the "basic model set" (i.e. set "c") of the Geneva isochrones 
1992). Temporal interpolation is a much 



(Schaller et al. 



more significant issue for post-MS isochrones than the PMS 
isochrones we have fitted in the past, as there are sharp dis- 
continuities in the rate of change of magnitude and colour 
with time, as exemplified by the MS turn-off. We there- 
fore use the code provided on the website to interpolate the 
isochrones to the appropriate age. We then convert from 
luminosity and effective temperature to colours and magni- 



tudes using the tables of Bessell et al. (1998), assuming the 



colours of Vega are zero (though V — 0.03). We also use 
Bessell et al's colour dependent extinction vectors. 

For some of the most luminous stars the gravities are 
rather low, and fall just outside the range of gravities given 



by Bessell et al. (1998). In these cases we extrapolate the 



models by simply setting the colour to that for the lowest 
available gravity. In these cases a linear extrapolation would 
be different by less than 0.001 mags, implying that the over- 
all error due to the extrapolation is much smaller than the 
uncertainties in colour. 

For reasons explained in Section [9~2| we used the Tycho- 
2 photometry for a Ori. In this case we have used the con- 
version given in Bessell ( 2000 ) to convert the Geneva-Bessell 
isochrones into the Tycho system. (H0g et al. 2000 state 



that the Tycho- 1 and Tycho- 2 systems should be identical.) 
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We used the reddening vector derived in M ayne &; Naylor] 
(2008). 



4 STATISTICS 



In Naylor &; Jeffries ( 2006 ) we introduced a solution to the 



long-standing problem of how to fit photometric data to 
isochronal models in colour-magnitude diagrams (CMDs). 
Whilst fitting an isochronal model (a curve) to a set of data 
points may at first appear to be a simple % 2 problem, the 
facts that the data points have uncertainties in two dimen- 
sions, and that the curve is smeared by binarity into a two- 
dimensional distribution, means a more sophisticated tech- 
nique is required. We have now used our solution to derive 



ages and distances for the young clusters NGC2547 (Nay- 
lor fe Jeffries|2006l) and NGC2169 fl Jeffries et al.|2007| ), and 
consistent distances to some of the best-studied star-forming 



regions ( [Mayne fe Nay lor 2008). In | Jeffries et al.| ( |20Q9| ) we 
derive distances to Vel OB2 and the association around 7 



Vel, and Joshi et al. ( 2008 ) use our technique for measuring 



the distance to NGC7419. 

Naylor & Jeffries! (120061) provide a rigorous development 
of r 2 , but for the purpose of understanding the improve- 
ments we have had to make to the method, a relatively 
simple intuitive interpretation gives a better insight into the 
problems. Figure [l] shows a typical fit of a dataset (shown as 
circled error bars) to a model (the colour scale). The model 
is a simulation of roughly a million stars (including binaries) 
using a specific age, metallicity, mass function and distance, 
which is then sampled onto a grid in colour-magnitude space. 
A fit in (for example) distance can be viewed as moving 
the model in the ^-direction until one obtains the strongest 
overlap between the model and the data. This overlap can 
be quantified for a single data point by taking the function 
which represents its position and uncertainties (normally a 
two-dimensional Gaussian), multiplying it on a gridpoint- 
by-gridpoint basis by the grid, and then summing the re- 
sulting values. If the grid is p(c, m) (where c and m are the 
colour and magnitude co-ordinates respectively) and the ith 
data point and its uncertainties Ui(c — Ci,m — mi) (where 
(ci,rrii) are its co-ordinates), then mathematically the over- 
lap is the integral of Uip over the entire space. The product 
of these integrals for all the data points will therefore reflect 
the overall overlap between the data points and the model, 
and so we define a statistic 
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whose minimum value corresponds to the best fit. 



In Naylor Sz Jeffries (2006) we showed that this defini- 
tion will, for models which are curves in (c, m) space, and 
only have uncertainties in the ra-axis, reduce to that for \ 2 • 
However, this is only the case if one chooses to multiply p by 
a normalisation factor which is dependent on the gradient 
of the isochrone. Unfortunately this normalisation factor be- 
comes infinite if the isochrone is vertical, and double-valued 
at any magnitude at which the isochrone is double valued. 
This means our x 2 -like normalisation will fail for the CMD 
fitting required here, because as one moves up the sequence 
towards bright magnitudes the isochrones become vertical, 
before finally switching to a negative gradient. Furthermore, 



if we wish to fit in U — B / B — V space, the isochrones are 
double- valued for certain values of B — V. 

In what follows, we therefore develop an alternative nor- 
malisation, which allows us to fit the data. In doing so we 
expose the limitations of an approximation we made when 
calculating the probability that the data are a good fit to 
the model. 



5 FITTING THE DATA 
5.1 The model CMD 

We must first create a probability density function to fit 



to the data. As in Naylor &; Jeffries (2006), we create this 



by simulating stars over the appropriate range of masses. 
For each star we choose a mass randomly from the Salpeter 
IMF, and if the star is a binary, we assign it a companion of 
a mass drawn from a uniform distribution between zero and 
the mass of the primary. The stellar model then provides a 
luminosity, gravity and effective temperature for each star, 
which we then convert into colour and magnitude using the 
appropriate bolometric corrections. If a binary companion is 
so low mass, or so cool, that it does not appear in the mod- 
els it is assigned a flux of zero. (Note that this assignment is 
a change from Naylor & Jeffries (2006), but has been used 



in all our subsequent work.) The value of each pixel in the 
image is then simply the number of stars whose colours and 
magnitudes lie within the pixel. We typically simulate 10 6 
stars, and for this work have used pixels of size 0.0025 mag- 
nitudes in each axis. This is half the value we have used in 
previous work, but is necessitated by the small uncertainties 
of the current data. We find the residual effects of the place- 
ment of pixel boundaries are much smaller (^0.005 mags in 
derived distance modulus) than the uncertainties in derived 
parameters. 



5.2 The normalisation of p 

Before proceeding further we must address the normalisation 



of the model image, p. In Naylor & Jeffries (2006) we used 
our x 2 -like normalisation which was a function of magni- 
tude. Here, we instead explore the results of a much simpler 
normalisation, setting the integral of p over the entire image 
to one. This raises the question of how faint magnitudes we 
must integrate down to. In fact, the strictly correct way to 
proceed would be to first multiply the image by the pho- 
tometric completeness function, such that below a certain 
magnitude p was zero, and then set the integral of what 
remains to one. Such a normalisation has an interesting, 
though subtle implication. When fitting for distance mod- 
ulus as the distance modulus increases, there is a decrease 
in the non-zero area of p, the region between the faintest 
observable absolute magnitude (as defined by the complete- 
ness function) and the brightest model star. Given that the 
integral over the model remains one, this means the value 
of any non-zero pixel will increase, implying that r 2 will de- 
crease, and hence the fit improve. This is actually the correct 
behaviour since it means that a model which fully populates 
the upper part of the sequence is better than one which does 
not. Practically, for our data, we can use a simpler normali- 
sation, where we make the integral between the faintest and 
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the brightest data points one. This means we have thrown 
away one possible source of information, but in practice this 
does not significantly affect the fits. 

Comparing the results obtained using this normalisa- 
with 



Naylor & Jeffries (2006) simply 

T" 



tion with that used in 

changes the values of r 2 in a given r z parameter grid by 
an additive factor; it does not change the best-fitting pa- 
rameters. This is at first a surprising result since we are 
changing the value of p in one part of the isochrone com- 
pared with another, which may appear as a weighting of the 
points. However, it should be remembered that adding the 
logarithms of the integrals in Equation [I] is equivalent to 
multiplying them together, so changing the relative values 
of p as a function of magnitude is a normalisation, not a 
weighting process. The only possibility for altering the best 
fit is if the length scale for changes in p is small compared 
to the size of an error bar. Then data points will drag the fit 
so that they lie in the higher-valued regions of p. Since the 
data point is more likely to originate in the higher- density 
part of the model this would again be the correct behaviour. 

Finally, it is important to note that we have no longer 
"normalised out" the mass function as we did in |Naylor fe| 
Jeffries (2006). Changing the mass function will change the 
value of r . In practice we have chosen to fix it such that 
dN/dM oc —2.35, which results in good fits to the models. 



5.3 The normalisation of U 

At the same time as considering the normalisation of the 
model, we should also consider that of the uncertainty func- 
tion (U in Equation [l]) . In % 2 fitting this is set such that 
the maximum value of U is always the same, so the high- 
est probability attainable is always the same, correspond- 
ing to a perfect fit, i.e. % 2 = 0. This is the normalisation 



we adopted in Naylor & Jeffries ( 2006| ). However, there is 
another obvious possibility, setting the integral of U to be 
one. This would have a very significant advantage in cases 
where the error bars seem to have been significantly under- 
estimated, and to obtain a good fit (i.e. a value of r 2 which 
corresponds to a Pr(r 2 ) of approximately 0.5) one has to 
add an extra uncertainty to U, in addition to those from 
the observations. This could well be due to mis-matches be- 
tween photometric systems. In such cases the procedure we 
have previously adopted has been to calculate r 2 , and then 
Pr(r 2 ) for increasing values of the added uncertainty, un- 
til Pr(r 2 ) exceeds 0.5. However, if one normalises U such 
that its integral is one, then conceptually one is comparing 
a model which includes the uncertainties with data points 
which are ^-function. One can, therefore, simply adjust the 
values of the uncertainties until one obtains the lowest value 
of r 2 . 

This normalisation has an additional conceptual advan- 
tage. In the case where the uncertainties are very small one 
can now approximate [/asa 2-dimensional ^-functions. This 
effectively removes the integral in Equation ^ and means 
one can evaluate r 2 by simply multiplying together the val- 
ues of p at the positions of the data points. 

We will refer to a normalisation where the integral of the 
models and the integrals of the uncertainty functions are all 
one as the natural normalisation. This cl early distinguishes 
it from the x 2_ like normalisation used in Naylor Jeffries 
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Figure 1. The data and best-fitting model for Cep OB3b. The 
colour scale is the model (p in Equation^) and the encircled error 
bars are the data. 



5.4 The fit 

Given we now have the correct normalisations we can now 
fit our example data, which is a sample from Cep OB3b de- 
scribed in detail in Section T9. 61 We calculated the extinction 
on a star- by-star basis as described in Section [5] and after 
correcting for it, searched in both age and distance modu- 
lus, evaluating Equation[l]at values of our fitting parameters 
which cover the range of interest. The resulting r 2 space is 
shown in Figure [2] The best fit, which lies at 10 Myr and 
a true distance modulus of 8.7 mags, is shown overlayed on 
the data in Figure [I] 

In some fits we find that there are data points which 
clearly do not lie on the sequence, and are presumably non- 
members. To deal with these objects we first fit the data 
with a variant of the "soft clipping" first described in Sec- 



tion 7.1 of Naylor & Jeffries ( 2006 ). We adapt this to the new 
normalisation by imposing a maximum r 2 for any one data 
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Figure 2. The r 2 grid for Cep OB3b. The contour is at the 68 
percent confidence level. 



point. The value used is the minimum value of r 2 amongst 
all the data points, plus a fixed value, normally 20. We im- 
plement this by calculating the probability corresponding to 
the imposed maximum r 2 and adding this to the calculated 
probabilities for each data point before calculating their r 2 
values. We then performed a second fit removing the data 
points which had r 2 values close to the clipping limit, with 
no clipping limit applied. 



6 IS THE MODEL A GOOD FIT? 

To test whether the model is a good fit, one must evaluate 
the chance of obtaining a given r 2 or below. One does this 
by calculating Pr(r 2 ), t he cumulative distributio n of the ex- 



pected value of r . In Naylor &; Jeffries (2006) we showed 



how to calculate this, for no free parameters, in such a way 
that it was insensitive to an incorrect choice of mass func- 
tion. We then suggested that one allow for free parameters 
by multiplying the r 2 axis of the distribution by (N — n)/N. 
Although our numerical simulations in Naylor & Jeffries 



0.6 



u 



0.2 



0.0 



\\ 
\\ 


\\ 
\\ 




\\ 


\\ 

V SOLID: x 2 - 


-LIKE 


\\ 
\\ 


\ DASHED: ) 


f (30 DoF) 


\\ 

SOLID: NATURAL V 


\\ 
\\ 
\\ 




DASHED: GAUSSIAN \ 

I V 


\ 

i Vi 





-50 



100 



150 



(2006) showed that the above approach to the free param- 



Figure 3. The probability of obtaining a given r 2 for a fit of 
30 data points to a main-sequence. This is the distribution one 
would obtain if one created a large number of datasets at a given 
distance modulus and extinction, and then "fitted" the data with 
the distance modulus and extinction fixed at their original values. 
The right-hand solid curve is for the x 2 -like normalisation, the 
left-hand solid curve for a natural normalisation where the sum 
of the probability over all colours is independent of magnitude. 
For comparison the dashed curves show the x 2 distribution for 30 
degrees of freedom and the Gaussian distribution for a = 60, with 
their expectation values shifted to match those for the x 2 -like and 
natural distributions respectively. 



eter problem may be approximately correct for the x 2_ like 
normalisation, it is straightforward to show that it cannot 
be correct in an arbitrary normalisation, such as the one 
described in Section |5.2| Consider a plot of the cumulative 
distribution of Pr(r 2 ) as a function of r 2 (Figure [3]). Chang- 
ing the normalisation of the model means multiplying p in 
Equation [l] by a constant. This has the effect of adding a 
constant to the values of r 2 as shown by the solid curves 
in Figure [3] Allowing for free parameters by scaling the r 2 
axis of the distribution by (N — n)/N would yield different 
values for the decrease in r 2 when adding extra parameters, 
depending on the normalisation. This cannot be correct, the 
decrease must be additive for the shape of the distribution 
to be invariant for a change in normalisation. 

There is an approximate solution to this problem, 
though, based on the fact that for CMD fitting, the distri- 
bution of Pr (t 2 ) is similar to Pr (x 2 ) 5 save a additive factor. 
The reason for this is that both distributions derive from 
the distribution of probability, and therefore r 2 , in the CMD 
plane. For a % 2 problem this distribution is a line, smeared 
by a one-dimensional Gaussian. For the r 2 CMD problem 
the distribution approximates to two sequences (those of 
single-stars and of equal-mass binaries) smeared by a two- 
dimensional Gaussian. It is, therefore, unsurprising that the 
resulting distributions of Pr (r 2 ) and Pr (% 2 ) are similar. So 
we could approximate Pr(r 2 ) by simply using the % 2 dis- 
tribution directly. However, for large values of N — n the 
differential form of the \ 2 distribution tends to a normal 
distribution whose mean is N — n and whose a is 2(N — n)Q 



It is interesting that the differential form of the r 2 distribution 
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This means we can allow for n free parameters by subtract- 
ing the expectation value from the distribution of t 2 , mul- 
tiplying the t 2 axis by (N — n)/N), and adding back the 
expectation value less n. We have implemented the latter 
approach, as it retains any asymmetry in the distribution. 
Applying this to the Cep Ob3b data results in a value of 
Pr(r 2 ) of 0.05. This is on the margins of acceptability, but 
no one datapoint is clearly discrepant. 



7 UNCERTAINTIES 

We have found a faster method for calculating the uncertain- 
ties than that presented in Naylor & Jeffries ( 2006 ) . The aim 
of the calculation is to place a contour in the r grid of Fig- 
ure [2] which represents a region within which the parameters 
lie with a given confidence. We can derive the uncertainties 
by first converting the values of r 2 in the grid into probabil- 
ity, and then integrating over the entire grid. We then divide 
this into the integral of the probabilities below progressively 
higher values of r 2 to obtain the cumulative r 2 distribution. 
We can then pick off values of r 2 at given confidence limits, 
and draw contours on the r 2 space. 

There are four practical issues which have to be solved 
when using this method. The first is that to carry out the 
integral one must multiply each pixel by its area. If the axes 
are linear then the the area of the pixels is the same, and the 
sum of the pixels will suffice, as we normalise by the integral 
over the whole area. However, if the age axis is logarithmic 
the simplest method is to multiply the probability by the 
age for that pixel, before performing the sum. 

The second problem is the underlying assumption that 
the model is correct. This means that the fitting to create the 
grid must be carried out using only those data points which 
are consistent with the model. So practically this means a 
second fit must be carried out excluding any points which 
the first fit clipped, without any further clipping (see Section 
5.4). Even so, this means one fit as opposed to fitting typ- 
ically 100 Monte-Carlo datasets for the previous technique, 
giving a speed improvement of a factor of 100. 

The third issue is that one must sum the grid out to in- 
finity. This is less demanding than it might at first appear. 
For example, if fitting a single data point in one dimension 
with Gaussian uncertainties, one only has to move ±3cr from 
the best fit to include 99.7 percent of the total probability, 
which is accurate enough for calculating a 95 percent confi- 
dence interval. Note, however, that the probability enclosed 
for a given a declines as the power of the number of observ- 
ables measured for each data point, and so if generalising r 2 
to many dimensions one would have to act with caution. 

The final issue is that the machine precision may be ex- 
hausted for some data points towards the edge of the r 2 grid, 
where the corresponding values of p are very close to zero. 
For example, if the smallest represent able number greater 
than zero is 1 x 10 -36 , the highest r 2 which can be obtained 



is approximately 168. Once the t 2 for any data point lies be- 
low this probability, the computer will calculate r 2 for the 
whole data set to be infinite. To flag such points in the grid 
we set them to a high value of r 2 (the number of data points 
times the r 2 resulting from a probability of twice the small- 
est representable number). Since such a r 2 is guaranteed to 
return a probability of zero, this works transparently in code 
which implements the new method of calculating confidence 
limits. However, we also note this number in the header of 
the grid file, so its meaning is clear if plots are made from the 
file. Applying this technique to the Cep OB3b data results 
in the contour shown in Figure [2] 

For the work here, the distance is a nuisance parameter, 
and we need to be able to quote an uncertainty in age alone. 
We therefore integrate the probability in Figure [2] over all 
distance modulii at each value of the age to create a run of 
probability with age. We then define a confidence limit as 
that region in age which integrates to give 68 percent of the 
probability, and which excludes equal integrals of probability 
above and below it. For Cep OB3b this gives an age range 
of 8.6 - 10.9 Myr. 



8 THE EXTINCTION 

Now having shown how we can fit for age, we must return 
to the question of the extinction. We follow an improved 
version of the two-strand approach developed in Mayne <S| 
Naylor| ([2008|) . We first attempt to fit the the U - B/B - V 



data with just the reddening as a free parameter. We can 
now (in contrast to |Mayne &; Naylor |2008| , test whether the 
model is a good description of the data. If it is, we assume 
the extinction is uniform, and apply the derived extinction 
to all the data points. If Pr(r 2 ) is too high, we conclude the 
extinction is non-uniform, and resort to deriving individual 
extinctions for each star by moving them along the (colour 
dependent) reddening vector until they reach the single-star 
U — B I B — V isochrone. This is essentially a modern version 
of the Q method of Johnson & Morgan (1953). As explained 
in |Mayne &; N aylor ( 2008 ) the disadvantage of this method 
is that it cannot allow for the fact the star may be a binary. 
This has the effect of narrowing the dereddened sequence 
in V/B — V space, hence our preference for the r 2 method 
where the extinction can be shown to be uniform. 



9 MAIN-SEQUENCE AGES 

We can now apply our technique to the rest of our sample 
of clusters and associations to derive MS ages. Each dataset 
we fit is given as an (electronic only) table as summarised in 
Table [TJ though we show the data for A Ori as an example 
in Table [2] 



(like the x 2 distribution) tends towards a Gaussian. This is a nat- 
ural consequence of the central limit theorem, since multiplying 
functions together and then taking the logarithm is equivalent 
to averaging their logarithms. Whilst the mean and width of the 
distribution are problem dependent, this may still provide a key 
to the solution in more general cases. 



9.1 NGC6530 



We used the data and uncertainties of Walker (1957), which 
are for a sample which is unbiased in colour and taken from 
an specific area of the cluster. To ensure we excluded the 
PMS, we selected only those stars blueward of B — F=0.28 
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Brown et al. (1994). The resulting r 2 contour does not close 




-0.4-0.3-0.2-0.1 
Intrinsic B-V 



0.3-0.2-0.1 
Intrinsic B-V 



Figure 4. The best-fitting model for NGC6530 with the age 
fixed at 2Myr (left) and the best fitting model with both age and 
distance as free parameters (right). Note how the five brightest 
datapoints are better fitted in the right-hand plot, and how the 
group of datapoints below them would have the be interpreted as 
binaries in the model on the left. The very faintest single stars 
also lie marginally closer to the single-star sequence in the better- 
fitting model. 



at low ages, and so we have only an upper limit on the age. 
We therefore quote (in Table [TJ the upper limit below which 
68 percent of the probability lies, but exclude this cluster 
from further analysis. 



9.3 NGC2264 



We used the photo-electric data of Walker (1956), as pre- 
sented in his Table 1. Fitting all stars blueward of B — V — 
for extinction in U — B / B — V space gives Pr(r 2 )=0.37, im- 
plying uniform extinction over the field. We then fitted in 
V IB - V and obtained a Pr(r 2 ) of 0.06. This is on the 
margins of acceptability, and there is a case that the two 
data points furthest redward from the sequence should be 
removed. However, in not doing so we simply enlarge our 
uncertainty estimate, and so are being conservative. 



9.4 A Ori (Collinder 69) 



We used the data from |Murdin &; Penston (1977) taking 
only those stars within half a degree of A Ori. We excluded 
objects with B — V > 0.2, which results in a sample which 
is almost complete blueward of this colour, and has no stars 
redward of (B — V)o = —0.04. After applying reddenings 
determined on a star- by-star basis, we obtained a value of 
Pr(r 2 ) of 0.52, provided we assumed the uncertainties were 



0.01 mags in V and 0.008 mags in B — V (Murdin & Penston 
1977, do not provide error bars) and removed two objects 
(HD36881 and HD36913) which appear to be non-members 
based on their position in the Vo/(B — V)o diagram. The 
resulting fit is shown on the right-hand side of Figure [5] 



and brighter than V = 13 and derived individual extinc- 
tions. The resulting fit is shown on the right-hand side of 
Figure [4] 



9.2 a Ori 



Our sample consisted of the members listed by [Sherry et al.| 
(2008} that are bl ueward of B — V = 0.1. We omitted the 
two stars noted by Sherry et al. (2008) as variable, and HD 



37333 which is above the MS, and probably a PMS star. 
We could not find a consistent Johnson UBV dataset for 
this cluster, and so used the Tycho-2 catalogue and its first 
supplement ( H0g et aL]|2000 ), although it does not contain 
a magnitude for a Ori C. In this dataset a combined mag- 
nitude is given for a Ori A and B. As a Ori A is itself a 
binary, we removed the effect of a Ori B on the combined 
magnitude by assuming the magnitude difference between 
the components is the mean of the values for the difference 
found from speckle ( Horch et al.||200T ) and adaptive optics 
( ten Brummelaar et al.||2000 ) work, a Ori B will have little 
effect on the combined colour. 

The disadvantage of using the Tycho-2 data is that we 
cannot determine the extinction, as there are no [/-band 
data. We therefore simply adopted E(B — V) = 0.06 from 



9.5 NGC2362 

We used the data and associated uncertainties for NGC2362 
from Johnson &; Morgan| ( |1953| ), which were taken as part of 
a programme to define what became the UBV system. We 
used only those stars blueward of B — V = 0.04 and excluded 



stars noted as non-members by Johnson & Morgan ( 1953 ) 



We also excluded the brightest star (r CMa) as it is clearly 
beyond the turnoff. Finally we found that star 36 gave a 
high r 2 in both U — B/B — V and V/B — V, and star 50 
in V/B — V, and so removed them from the fit as well. We 
then measured a global extinction from the U — B / B — V 
diagram, before fitting in V/B — V. 



9.6 Cep OB3b 



Blaauw et al. ( 1959 ) carried out a photometric survey of 



stars of spectral type AO and earlier identified from objec- 



tive prism plates. We take the membership list from |Pozzo| 
(2001), but exclude BHJ11 for which the measurement is a 
combined light measurement for a rather wide Am = 2.5 
binary. We de-reddened this sample on a star- by- star basis. 
Using the uncertainties quoted in the paper, we obtain a 
just about acceptable value of Pr(r 2 )=0.05. 
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Figure 5. The best-fitting model for A Ori with the age fixed at 
3Myr (left) and the best fitting model with both age and distance 
as free parameters (right). The brightest three stars are clearly 
better fitted by the model on the right. 
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Figure 6. The best-fitting model for the ONC with the age 
fixed at PMS age of 3Myr (left) and the best fitting model with 
both age and distance as free parameters (right). The 68 percent 
confidence interval for the best-fitting age just encompasses 3Myr, 
and so we expect the improvement in the fit from left to right to 
be only marginal. We can see the improvement in the fit is due 
entirely to the improvement in the fit for the brightest star, and 
even that is at the expense of a worse fit for the second brightest 
star. Thus the conclusion that the statistics drive us to, that the 
improvement is marginal, seems reasonable. 



9.7 The Environs of the Orion Nebula Cluster 



9.8 NGC2547 



Our data and uncertainties are taken from Walker (1969), 
who aimed to obtain photometry for as many stars as possi- 
ble within the outline of the dark cloud, since this area will 
be the least contaminated by background stars. We removed 
stars redward of B — V=0.0, those marked as variables or 
visual doubles, and three stars which lie away from the se- 
quence in the U — B / B — V diagram. We fitted the data for 
a single extinction in U — B / B — V space, and after remov- 
ing two outliers in r 2 obtained a good fit with Pr(r 2 )=0.46. 
At first this may sound counter-intuitive, since it implies 
uniform extinction, yet it is well known that the extinction 
of ONC members is highly variable. In fact it seems this 
only applies to stars in the central cluster. We then fitted in 
V / B — V to obtain the results in Table [I The resulting fit 
is shown on the right-hand side of Figure 6] 



We took the photometry from Claria| ([1982 ) , and used an 
uncertainty of 0.02 mags in both magnitude and colours. 
This is an estimate for U — B and corresponds to the devi- 



ations for single observations derived by Claria (1982) from 
comparison with the data oflFernie (1959). Although Claria 



(119821) often has four observations per star we prefer to take 



the view that the uncertainty represents the difference be- 
tween the photometric systems. We use the membership list 



of Claria ( 1982 ), which is based on proper motions, photom- 
etry and spectroscopy, and select only stars with B — V < 0.1 
to ensure we exclude PMS stars. We found that if we in- 
cluded star 40, which appears to sit just below the MS we 
obtained an unacceptably low value of Pr(r 2 ). Furthermore 
this star is right on the edge of the proper motion distribu- 
tion of the bulk of the members, so we excluded it. 
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9.9 IC2602 



We used the data of Eggen ( 1972 ), excluding stars with B 



V > 0.0, and HD93163, which lies away from the sequence. 
Unfortunately Eggen (1972) does not provide uncertainties, 
but we found a single extinction would yield Pr(r 2 )=0.79 for 
uncertainties of 0.014 and 0.014 mags in B — V and U — B 
respectively, which suggests the extinction is uniform. Using 
that extinction, and an uncertainty in V of 0.025 mags gives 
Pr(r 2 )=0.27 when fitted inVvsB — V. 

9.10 The Pleiades 



We again used the data and memberships from Johns on~fe| 
Morgan| ( |1953| ). The U - B/B -V diagram, especially the 
region where the gradient is reversed, shows that there is 
variable extinction to this cluster. We therefore dereddened 
the data on a star- by-star basis, which limits us to B — V < 
0.0. Before fitting we also excluded Hertzsprung 371 (which 
appears to be reddened). 



10 PRE-MAIN-SEQUENCE AGES 

For the PMS ages we require a set of consistent ages, and 
we therefore adopt the ages of |Mayne &; Naylor| ( |2008| , with 
the following exceptions. 

10.1 NGC2547, IC2602 and Cep OB3b 



We t ake th e PMS age for NGC2547 from |Naylor JeT] 
fries (2006) as 38.5Myr, which is derived from isochrone fit- 



ting, though also agrees with the Lithium depletion age. The 



age for IC2602 (25Myr) is taken from |Stauffer et al.| fll997| ), 
which again is based on isochrone fitting to the PMS stars. 
We use a PMS age of 4.5Myr for Cep OB3b, which is from 



Littlefair at al (in prep) , but is based on the system of Mayne 
& Naylor|([2008). 



10.2 The Environs of the Orion Nebula Cluster 

The position on the sky of our MS sample is shown in Fig- 
ure |7| along with the positions of the sample of Hillenbrand 
( |1997| , which represents stars in the ONC itself, and the 
flanking fields of | Ramirez et aL (2004). Given the distribu- 
tion of stars, it is clear that the PMS age we should use is 



that of the flanking fields. Although Ramirez et al. (2004) 



calculate this, they do so on the assumption that the ONC is 



470pc away, whilst a more modern estimate is 400pc (Mayne 
fe Naylor||2008l and references therein). In the V /(V - I)o 



diagram Ramirez et al. ( 2004| place the flanking fields 0.3 
mags above NGC2264. Correcting the distance to 400pc will 
bring the flanking fields PMS to the same magnitude as that 
of NGC2264, and therefore to an age of 3Myr on the scale 



of Mayne fe Naylor (2008) 
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Figure 7. The positions of stars in the vicinity of the Orion 
Nebula Cluster. The dots in the central region are one in five of the 
stars from Hillendbrand (1997). The dots around the periphery 
are X-ray sources from Ramirez et al. (2004). The filled circles 
are the main-sequence sample. 



Table 2. An sample of Tables 2-11, the fitted dataset for A Ori. 
Table [l] gives the number of the electronic table for each dataset, 
along with the reference for the star numbering system. As shown 
here, for each cluster we give, along with the uncertainties, the 
fitted V, B — V and U — B. In the case of groups where extinctions 
were derived on a star-by-star basis these are the reddening and 
extinction free values, and the E(B — V) used is given in the last 
column. 



Star 
number 

36822 
36861 
36862 
36894 
36895 
37034 
37035 
37051 
37110 
245140 
245168 
245185 
245203 



V 



mag 

4.036 
3.346 
4.710 
8.662 
6.512 
8.976 
8.249 
8.699 
8.657 
8.568 
9.077 
9.313 
6.972 



0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 



B - 
mag 

-0.270 
-0.271 
-0.241 
-0.086 
-0.188 
-0.069 
-0.141 
-0.074 
-0.097 
-0.107 
-0.057 
-0.050 
-0.188 



• V 



0.008 
0.008 
0.008 
0.008 
0.008 
0.008 
0.008 
0.008 
0.008 
0.008 
0.008 
0.008 
0.008 



u - 

mag 

-1.047 
-1.049 
-0.953 
-0.295 
-0.745 
-0.216 
-0.562 
-0.239 
-0.346 
-0.398 
-0.173 
-0.153 
-0.745 



0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 



E(B-V) 



0.120 
0.061 
0.281 
0.036 
0.068 
0.109 
0.121 
0.114 
0.097 
0.217 
0.177 
0.190 
0.158 



we emphasise that those of Mayne & Naylor (2008 ) are to be 
preferred. We plot PMS against MS age in Figure [8] Whilst 
the PMS and MS ages for individual clusters may agree to 
within the uncertainties, the average of the MS ages is sig- 
nificantly older than the average of the PMS ages. If we 
take only those clusters and associations less than lOMyr 
old, the MS ages are, on average, a factor two larger. The 
issue is clearly which of these age scales is correct. 



11 DISCUSSION 

We collect together our measurements of the ages of the 
groups and clusters in Table [I] along with the other pa- 
rameters from our fits. For completeness we include the dis- 
tances, though as these are derived from two-parameter fits 



11.1 Are the main-sequence ages incorrect? 

Explanations as to why the MS ages may be incorrect fall 
into two groups, those associated with the statistical tech- 
niques and those associated with the models. We can rule 
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Table 1. Main-sequence and pre-main-sequence ages. 



Cluster or 


PMS Age 


MS Age (Myr) 


Pr(r 2 ) 


Distance 


E(B-V) 


Data 


Group 


(Myr) 


Best fit 


68% confidence 




Modulus 




table 


A Ori 


3 


6.6 


5.8-7.5 


0.52 


8.05 


0.12 1 


2 


NGC6530 


2 


5.5 


4.9-6.1 


0.84 


10.76 


0.33 1 


3 


NGC2264 


3 


5.5 


2.4-6.0 


0.06 


9.24 


0.05 


4 


ONC 


3 


5.0 


2.8-5.2 


0.10 


7.92 


0.03 


5 


a Ori 


3 


0.4 


< 6.6 


0.26 


7.98 


0.06 2 


6 


NGC2362 


4.5 


9.1 


5.4-12 


0.08 


10.71 


0.08 


7 


CepOB3b 


4.5 


10 


8.6-10.9 


0.05 


8.72 


0.89 1 


8 


IC2602 


25 


44 


28-62 


0.27 


5.88 


0.02 


9 


NGC2547 


38 


48 


27-62 


0.13 


8.03 


0.04 


10 


Pleiades 




115 


104-117 


0.81 


5.35 


0.02 1 


11 



1 Median from individual extinctions. 
2 From |Brown et al.||l994| . 



Numbering 
system reference 



HP catalogue 
" 1957) 



Walker 
Walker |T956> 
^unfll^ 



Johnson 


|l950 


) 


Blaauw et al. 


1959 



|Claria| jl982 




|Hertzsprung 


(1947 




PMS CONTRACTION AGE (MYR) PMS CONTRACTION AGE (MYR) 



Figure 8. The main-sequence and pre-main-sequence ages for 
our sample. The groups at PMS ages of 3 and 4.5Myr have been 
separated slightly in age to aid visibility. 



Figure 9. As ngure[8]but with pre-main-sequence ages calculated 
without the brightest star for each dataset. 



out problems with the fitting procedure by comparing our 



ages with those obtained by Meynet et al. ( 1993| . They use 
similar isochrones to the ones presented here and measure 
the age of the Pleiades as 100 Myr and the environs of the 
ONC as 4Myr. Both these ages are compatible with those we 
measure, suggesting our technique gives similar ages to "by 
eye" fitting. Equally importantly, we match the lithium de- 



pletion age for the NGC2547 (34-36Myr | Jeffries fe Oliveira 
2005 ) and are very close to the depletion age for the Pleiades 
(125-130Myr Stauffer et alTl998) . 

To check that our uncertainties are at least reasonable 
we tested how the result changes if one star is removed from 
each fit. As one might expect, the brightest star in the fit 
provides the tightest limits on the age. We therefore removed 
the brightest star from each dataset and replotted Figure [8] 
As Figure [9] shows, the result remains clear, though as one 
might expect the error bars are larger, and one more dataset 
(the ONC) returns an upper limit for the age. This strongly 
suggests that our uncertainty estimates are reasonable, and 
the result is robust. This experiment also shows what the 
effect might be of a non-member being included in the fit. 
Were a non-member very far from the fitted sequence it 
would have been clipped out by the procedure described at 



the end of Section |5.4[ Were it close to the sequence, then 
it could deviate the fit sufficiently to have a reasonable r 2 , 
but then would only change the best fit by a small amount, 
similar to the effect of removing a data point. 

As a final check of the uncertainties, in Figures^] to [6] we 
plot the data over the best fitting models if the age is fixed 
at the PMS age (left) or left as a free parameter (right). The 
first two examples (NGC6530 and A Ori) are ones where the 
PMS age lies far outside the 68 percent confidence region for 
the upper-main-sequence age. As one expects, we see that 
the brightest stars lie to the right of the model when the age 
is fixed at the PMS age. Our final example, the ONC, is one 
where the PMS age lies almost exactly on the edge of the 68 
percent confidence limit. Here the improvement in the fit is, 
as it should be, marginal. Although such comparisons with 
our expectations are at best subjective, that they fit with 
our expectations adds to our confidence in the result. When 
combined with the experiment of missing out the brightest 
datapoint, we have a strong case that our uncertainties are 
correct, and the result is robust. 

The obvious problems with the models are the absence 
of rotation, uncertainties as to the mass- loss rates, and the 
treatment of convective core overshoot. Figure 9 of Meynet 
& Maeder (2000) shows that if the stars were rotating, and 
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we fitted them with isochrones for stationary stars, the re- 
sulting ages would be too young by about 10 percent. This 
therefore exacerbates the discrepancy between the PMS and 
MS ages. 

All modern models include a degree of core overshoot, 
which has the effect of mixing more hydrogen into the core, 
and hence lengthening the MS lifetime. Naively, models with 
no overshoot will have shorter MS lifetimes than those used 
here, by roughly the decrease in available hydrogen (per- 
haps 20-40 percent), which is of the right order to bring the 
MS and PMS ages back into agreement. However, our CMD 
fitting does not measure lifetime on the MS, but how far 
from the ZAMS a star at a given luminosity (not mass) has 
moved. A close comparison of Figures 4 and 5 of Maeder & 
Mermilliod| ( |198l| shows that for the youngest ages they cal- 
culate (25Myr) the difference in the position of the isochrone 
corresponds to an age difference of around 5 percent. 

Mass-loss rates for early-type stars are uncertain, and 
so in addition to using the Geneva models with the standard 
mass- loss rates (set "c"of Mae der &; Meynet||1994 ) we also 
tried the higher mass-loss rate, set "e". Comparison of the 
resulting isochrones for the masses and ages we are inter- 
ested in shows differences in colour which are too small to 
affect our results. 

Finally, we have tested the effect of using different MS 
models. As an alternative to the Geneva models with the 



Bessell et al. (1998} conversions to colour and magnitude we 



used the conversions presented with the isochrones in |Leje-| 
une &; Schaerer (2001). We obtained ages somewhat older 



than those from the Geneva-Bessell models, exacerbating 
the age difference problem. More importantly, the values of 
Pr(r 2 ) are much worse than those for the Geneva-Bessell 
models, typically around 0.01 or 0.001, showing that these 
models can be ruled out as good descriptions of the data. 
To test whether this is the interior models or the atmo- 
spheres, we fitted the data to the Padova models (IGirardil 



et al.||2002 ) but with the same model atmospheres ( |Bessell| 
et al. 111998) as we used for the Geneva-Bessell models. We 



find this gives slightly younger ages (a factor 1.5 older than 
the PMS ages in the range l-10Myr), but very similar values 
of Pr(r 2 ) to the Geneva-Bessell models. In summary, there- 
fore, our fitting gives strong support for the Bes sell et al.| 
(1998) conversions, and there is only a weak effect from the 



els. Our conclusion, therefore, is that the MS age scale is 
probably the correct one. 



11.3 Implications of lengthening the PMS 
timescale 

Before discussing the implications of a longer timescale, we 
should be wary of over interpreting Figure [8] Whilst it 
clearly shows a discrepancy between mean PMS and MS 
ages, the error bars for individual data points are large. All 
we can say with any certainty is that there is a difference 
of approximately a factor two at PMS ages of 3Myr. By 
30Myr our data are consistent with the age scales matching, 
though a difference of a factor 1.5 is still, in the statisti- 
cal sense, likely. We therefore limit ourselves to discussing 
the implications of a lengthening of the timescales in the 
l-10Myr PMS age range. Even here, however, we find there 
are problems it might solve. 

There is a long-standing issue that the observed 
timescale for the dissipation of proto-stellar discs (3Myr; 
|Haisch et al.|2 001 ) may be shorter than the time required by 
the models for planet formation (lOMyr; Pollack et al. 1996). 
In recent years there has been significant effort to find mech- 
anisms which will shorten the planet forming timescales. 
Whilst a case can be made that this problem has been solved 
( Mordasi ni et al.|2Q 08 ) , there is a view that significant prob- 
lems remain (see, for example, the introductory sections 
of |Ayliffe fe Bate||2009[ |Dodson-Robinson et al.||2008| ). A 
fair summary is probably that whilst there are mechanisms 
which could shorten the timescale, such as dus t settling (|Hu- 
|bickyj et al. 2005) and planetary migration (Alibert et al. 
|2005| ), the uncertainties in the physics remain such that it is 
not clear they do. Our result offers an interesting alternative 
solution. If the clusters used to measure the disc dissipation 
timescale are 50-100 percent older than previously thought, 



there may be no contradiction with the Pollack et al. ( 1996) 
timescale. 



interior models, which can explain some, but not all of, the 
age discrepancy. 



Jeffr ies et al.| ( |2007 ) point out that there is a lack of 
clusters in the age range 5-30Myr. Revising the age scale 
in the way suggested by the MS fitting would move clusters 
from the youngest ages into this age range. Furthermore if 
the age scales come back into register at around 30Myr, as 
Figure [8] suggests they might, there would not be a compen- 
sating movement out of the 5-30Myr range, leading to an 
increased number of clusters at these ages. 



11.2 Are the pre-main-sequence ages incorrect? 

The PMS ages are much less robust than the MS ones. We 
have adopted the PMS age scale of [Mayne & Naylor| (|2008 ) . 



However, as Mayne &; Naylor (2008) and Mayne et al. (2007) 



make clear, the primary aim of this scale is an age ordering. 
The age scale itself is rather arbitrary, though was chosen 
to match as closely as possible the commonly quoted ages 
for the young groups. The problem is that there is no single 
PMS age scale, a point nicely illustrated in Jeff ries et al.| 
(2009). They show that the 7 Vel association could have a 
PMS age between 5 and 15Myr depending on which PMS 
models are used, and which part of the sequence is consid- 
ered. They estimate that the association is about 7Myr old 



on the Mayne & Naylor (2008) scale, so doubling the ages of 



these young associations is consistent with some PMS mod- 



12 CONCLUSIONS 

We have shown that there is a systematic difference between 
the ages of clusters and associations measured from the MS 
and ages commonly used which are based on the PMS. The 
difference is in the sense that the MS ages are a factor 1.5-2.0 
greater than the PMS ages in the age range 2-5Myr (on the 
PMS scale). The most straightforward solution is to adopt 
the MS age scale, as there are PMS models which fit with 
the longer timescale. Adopting the longer timescale offers a 
solution to the problem that the lifetimes of discs around 
stars (3-5Myr on the PMS age scale) are shorter than the 
time taken to form planets, and to the apparent absence of 
clusters in the 5-30Myr age range. 

Finally we should be clear that although we favour the 
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age-scale given by MS fitting, we are not recommending it 
as a method for deriving ages for individual clusters and 
associations. As Figure [8] and Table [I] make clear, the uncer- 
tainties for individual groups are large. Nor can we at this 
point make any clear recommendation as how one should 
reflect this result when quoting PMS ages. Whilst it is clear 
that the youngest ages need to be increased, how far down 
the age scale that should be propagated is unclear. We there- 
fore continue to commend the Mayne et al/Mayne & Nay- 
lor age ordering, though recommend that if these ages are 
quoted one states clearly that they are on the Mayne et 
al/Mayne & Naylor scale. If absolute ages are required for 
clusters younger than lOMyr for comparison with other data 
we recommend multiplying the Mayne et al/Mayne & Nay- 
lor values by 1.5 and quoting the age scale as originating 
from this paper. 
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